Entanglement and dephasing of quantum dissipative systems 
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The von Neumann entropy of various quantum dissipative models is calculated in order to dis- 
cuss the entanglement properties of these systems. First, integrable quantum dissipative models 
are discussed, i.e., the quantum Brownian motion and the quantum harmònic oscillator. In case 
of the free particle, the related entanglement of formation shows no non-analyticity. In case of 
the dissipative harmònic oscillator, there is a non-analyticity at the transition of underdamped to 
overdamped oscillations. We argue that this might be a general property of dissipative systems. We 
show that similar features arise in the dissipative two level system and study different regimes using 
sub-Ohmic, Ohmic and and super-Ohmic baths, within a scaling approach. 



I. INTRODUCTION 

The operation of a quantum computer requires a care- 
ful control of the interaction between the system and its 
macroscopic environment. The resulting entanglement 
between the system's degrees of freedom and the reser- 
voir has been a recurrent tòpic since the formulation of 
quantum mechanics, as it is relevant to the analysis of 
the measurement process 0, The loss of coherence 
due to the interaction of a quantum system and its en- 
vironment was also extensively studied due to its impor- 
tance to Macroscopic Quantum Tunneling and related ef- 
fects 0. Theoretical research on Macroscopic Quantum 
Tunneling lcad, among other rcsults, to the formulation 
of a canonical model for the analysis of a quantum sys- 
tem interacting with a macroscopic environment, the so 
called Caldeira-Leggett model Q. It can be shown that 
this harmònic model (see below) describes correctly the 
low energy features of a system which, in the classical 
limit, undergoes ohmic dissipation (linear friction). It 
can be extended to systems with more complicated, non 
linear, dissipative properties (usually called sub-Ohmic 
and super-Ohmic, see below) 

In relation to the ongoing research on entanglement, a 
recent interesting development is the analysis of the en- 
hancement of entanglement in a system near a quantum 
critical point 0, @- The original systems under study 
were the transverse Ising model and the XY model, but 
also other models which exhibit a quantum phase transi- 
tion were later investigated in this direction, as e.g. the 
Lipkin-Meshkov-Glick model [9l[ïo|. 

A connection between previous research on Macro- 
scopic Quantum Tunneling and entanglement near quan- 
tum critical points is starting to emerge EU. It is in- 
teresting to remark that some of the simplest systems 
which show a non trivial quantum critical point is the 
dissipative two level system [5( and related models, like 
the Kondo model [T^ . Similar models describe quantum 
fluctuations in Josephson junctions |l3j or tunneling be- 
tween Luttinger líquids 14] . It is alrcady known that, 
even for the ground state of simple models like the dissi- 
pative harmònic oscillator non-trivial entanglement prop- 
erties can be expected as was already commented on in 



Ref. Entanglement energètics at zero temperature 

was investigated in Ref. 0|. A number of properties 
of the entanglement in the Caldeira-Leggett model and 
related models remain, however, unexplored. 

The models studied here describe a quantum system 
characterized by a small number of degrees of freedom 
coupled to a macroscopic reservoir. These models show 
a crossover between different regimes, or even exhibit a 
quantum critical point. As this behavior is induced by 
the presence of a reservoir with a large number of de- 
grees of freedom, they can also be considered as a model 
of dephasing and loss of quantum coherence. It is worth 
noting that there is a close connection between models 
describing impurities coupled to a reservoir, and strongly 
correlated systems near a quantum critical point, as ev- 
idenced by Dynamical Mean Field Theory ^J- I n the 
limit of large coordination, the properties of an homoge- 
neous system can be reduced to those of an impurity in- 
teracting with an appropriately chosen reservoir. Hence, 
in the limit of large coordination the entanglement be- 
tween the quantum system and the reservoir near a phase 
transition can be mapped onto the entanglement which 
develops in an homogeneous system near a quantum crit- 
ical point. 

The measure of entanglement used in the original pa- 
pers is the the concurrence introduced by Wooters [l8| 
Alternatively, the von Neuman entro py o f macroscopic 
(contiguous) subsystems can be used [Ïq. A non-local 
measure of entanglement was employed in the study of 
the Affleck-Kennedy-Lieb-Tasaki (AKLT) model pol·luí]. 

Our previous work |ll| showed that the main non- 
analyticity of the concurrence arises at the transition of 
coherent to incoherent tunneling. At the actual quan- 
tum phase transition, we only found a much weaker 
non-analyticity associated with the existence of the 
Kosterlitz-Thouless weakly non analytical features. In 
the same way, the phase transition of the transverse Ising 
model discussed in Refs. 0, @ can be interpreted as a 
transition where coherence is lost due to the emergence of 
a localized state at the transition. We thus assume that 
the loss of coherence might be more important to see 
non-analyticities in the entanglement of a system than 
the actual phase transition. 
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In the first part of this article, we test our assump- 
tion using two integrable quantum dissipative models, 
the dissipative free particle - that is, the Caldeira-Leggett 
model - and the dissipative quantum harmònic oscillator. 
These models do not exhibit a quantum phase transition, 
but, in the latter case there is a transition from under- 
damped to overdamped oscillations at some critical cou- 
pling strength. As measure of entanglement we use the 
von Neuman entropy of the subsystem, defined using the 
reduced density matrix, p Al obtained by tracing out the 
bath degrees of freedom of the ground state: 

Ety) = -Tr(p A \np A ) , p A = Tr B {\Í>)(il>\) (1) 

In the second part of the paper, we make the same 
analysis for the spin-boson model on the basis of a scaling 
approach for the free energy. For super-Ohmic baths, the 
model shows no phase transition whereas for Ohmic and 
sub-Ohmic baths, there is a transition from localized to 
non-localized behavior. Again, we focus the discussion 
on the transition from coherent to incoherent oscillation 
which exists for Ohmic dissipation, but is also present for 
certain non-Ohmic environmcnts. 

It is finally worth noting that a more mathematical 
analysis of some problems related to the entanglement in 
the dissipative harmònic oscillator can be found in Refs. 
[22I I2ÜI ]. Sub-Ohmic environments may be relevant to 
the description no Gaussian effects in qubits coupled to 
external environments, see Í24L Í25Í . 



Leggett and further investigated by Hakim and Ambe- 
gaokar [2(| |27| • The latter authors obtained the reduced 
density matrix via diagonalization of the Hamiltonian: 



(x\ PA \x') =e- a{x - x '^/L 
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(5) 



where rj denotes the friction coefhcient and u> c is the cut- 
off frequency of the bath. Furthermore, L — > 00 denotes 
the system size and in contrary to the use of Eq. [SJin Ref. 
|27| , here the normalization is crucial to assure Trp A = 1 . 

In order to calculate the entropy of the system, we 
Taylor expand the logarithm: 



x ^ (1~paY 
\np A = - 2^ 



71=1 



Further we have 



(x\p A \x') = 
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proved by induction. With the idcntity 



(0) 



(7) 



(8) 



we thus obtain for the specific entropy (for general di- 
mension d) 



II. EXACTLY SOLVABLE DISSIPATIVE 
SYSTEMS 

Modcling the environment by a set of harmònic oscil- 
lators the general integrable model is described by 
the following Hamiltonian: 



„ P 2 . w o 2 . ÍPl , 1 2 / \ 



(2) 



S= -(ln(ai 2 ) + l-ln7r). 



(9) 



Comparing the above result with the entropy of a particle 
in a canonical ensemble, we identify a ~ A -2 cx T with 
A denoting the thermal de Broglie wavelength and T the 
temperature of the canonical ensemble. 

Notice that the entropy of a free dissipative particle 
shows no non-analyticity. 



The operators obey the canonical commutation relations 
which read (fi = 1) 



(3) 



The coupling of the system to the bath is completely 
determined by the spectral function 



(4) 



In the following, we will consider a Ohmic bath with 
J(u>) = i]U) for ui <C u> c and J(u>) — for u uj c , u c 
being the cutoff frequency. 



A. Caldeira-Leggett model 

Let us first consider the free dissipative particle, i.e., we 
set lüo = 0. The model was introduced by Caldeira and 



B. Dissipative harmònic oscillator 

We now include the harmònic potential, i.e., loq =/= 0. 
The reduced densitymatrix of the damped harmònic os- 
cillator is given by |(| 



{x\p A \x 




The above expression is deduced such that the correct 
variances for position and momcntum are obtained. At 
T = the expectation vàlues are given by 



<? 2 > - ^-/(«) 
(p 2 )=c 2 (l-2 K 2 )< (? 2 )-,- 



_ ,2, -i .,..2; 1 2w 0« j n / kV 

^0, 



(11) 

(12) 



with k — t]/2ujq and 



l ln [( K + V^Ï)/( K -V^Ï)] 
/(«) = ; j = 1 ■ ( 13 ) 

The parameter k represents the friction parameter and 
the system experiences a crossover from coherent to in- 
coherent oscillations for k = 1. 

Taylor expanding the logarithm of the entropy, Eq. 
©, we need to evaluate the general n-dimensional inte- 
gral 



dx\..dx n exp 



r n/2 



VdetA 



(14) 



where A is given by the translationally invariant tight- 
binding matrix with = 2(a + 6), Ai + u — = 
— (a — b) (n+1 = 1) and zero otherwise. The determinant 
of the matrix is given by its eigenvalues and reads 



detA = (2a) n (l - b/a) n J] 



1 



2b 



cos k„ 



(15) 
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FIG. 1: (Color online). The entropy S of the dissipative os- 
cillator (full) and the dissipative two-level system (dashed) 
with Ohmic coupling as function of the dimensionless cou- 
pling strength a. Inset: The derivative of the entropy S' with 
respect to the coupling strength a is shown as function of a. 



with k m = 2?:m/n. Considering the n-dimensional trans- 
lationally invariant, but non-hermitian matrix Aíj — 1, 

= 1 — e (n+1 = 1) and zero otherwise, one obtains 
the following formula: 



n 



i 



2(1 



COS kn 



n\2 



(1 -(! -£)") 
2 n {\ 



(16) 



For lüc/loq 3> 1, we have 
a/b = 4(q 2 )(p 2 ) 



/(*) 



4k 



(1-2^)/(k) + — ln 

7T \ CJ 



> 1. (17) 



In this limit, we can thus set £ 2 = 4&/a <gC 1 and the 
?i-dimensional integral can be approximated to yield 



dx{x\p 7 ) L \x) 



1 - (1 - e) n 



(18) 



with i = e\f\ — ej y/\ — e 2 /4. Expanding the denomina- 
tor as geometrical series, we have for the entropy 



S 



-lne+ 4 ln (! 

e e 1 



(19) 



In the limit e ss e <C 1, the leading behavior of the 
entropy is given by S ~ ln(a/ò). We thus find a non- 
analyticity at k = 1, the point of the crossover of incoher- 
ent to coherent oscillations. The leading behavior of Eq. 
(I19JI is plotted in Fig. \ï\a,s function of the dimensionless 
coupling strength a = (7^77/(2^) with the characteristic 
length scale q — l/y/UJo for üj c /üjo — 100 (full line). In 
the inset, the non-analyticity of the derivative of the en- 
tropy with respect to the coupling strength S' = d a S at 
a = X/t: can be seen. 



III. SPIN-BOSON MODEL 

A prominent dissipative model is given by the spin- 
boson model or dissipative two-level system (TLS). The 
Hamiltonian without bias reads 



H = ^-a x + u k b\b k + <T z J2^(b k + b{ 



At 



(20) 



The operators b^ resemble the bath degrees of freedom 
and a x , a y , a z denote the Pauli spin matrices. They 
obey the canonical commutation relations and the spin- 
1/2 àlgebra, respectively. 

The coupling constants A^ are parameterized by the 
spectral function 



J{íü) 



Afe«5(w - uj k ). 



(21) 



In the relevant low-energy regime, the spectral function 
is generally parameterized as a power-law, i.e., J(u>) oc 
2etLú s K.Q~ s where a denotes the coupling constant, s the 
bath type and Ao the cutoff-frequency. 

The general reduced density matrix of the spin-boson 
model is given by 



PA 



2 



1-((T Z 



(22) 



Since there is no symmetry breaking field in the above 
Hamiltonian, we have (a z ) — 0. The eigenvalues are thus 
given by A± = (1 ± (cr x ))/2 and the entropy reads 



S=- 



ln((l-(a x ) 2 ) /4) + (01n(i±^ 



(23) 
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The value of (a x ), at zero temperature, is given by 

/ i n dE 
^ = 2 dA 



(24) 



where E is the energy of the ground-state. To obtain the 
ground-state energy, a scaling analysis for the free energy 
at arbitrary temperature is considered, see the appendix. 
In the following, we use this approach to calculate E(A ) 
and (a x ) which will set the basis of our discussion on the 
entanglement properties of the spin-boson model. 



A. Ohmic dissipation 

In the Ohmic case (s = 1), there is a phase transi- 
tion at zero temperature at the critical coupling strength 
a = 1 [2^, [23. The transition is also manifested 
in the renormalized tunnel element A ren , i.e., A rcn = 
Ao(A /A ) q/(1 ^ q) for a < 1 and A rcn = for a > 1. 

The free energy is determined by (see the appendix) 



F = 



cm 



dA . 



(25) 



The ground state energy E can then be written as 



E = < 



Al » /AoV 
A7 ~ A ° V^J 



Ao 



Ao" 



< a < \ 

OL=\ 

\ < a < 1 
a > 1 



(26) 



where C is a numerical constant. For a 
have 



dln((7 x ) 



eia 

For a = 1 — e, we have 
dln(a x ) 



1 

oc -. 

a=l/2±e e 



oc ln — 

Q=l-e \ Aq 



1/2 ± e, we 
(27) 

(28) 



The non-analyticity around a — 1 is thus far weaker than 
around a — 1/2. This non-analyticity is also present in 
the entropy as can be seen from the expression Eq. (|22|l . 

The entropy S of the dissipative two-level system with 
Ohmic coupling is plotted in Fig. ^as function of the di- 
mensionless coupling strength a for w c / Ao = 100 (dashed 
line). The inset shows the derivative of the entropy with 
respect to the coupling strength. The entropy quickly 
saturates after the transition from coherent to incoher- 
ent oscillations at a = 1/2, but the non-analyticity of 
Eq. I|27|l cannot be seen on this seale. 



B. Non-Ohmic dissipation 

) The calculation of E(Aq) and (a x ) can be extended 
to the spin-boson model with non-Ohmic dissipation (s 7^ 
1). In general, the dependence of the effective tunneling 
term on the cutoff, A(A) is: 



A(A) = A exp 





f3H 







(29) 



with the spectral function given in Eq. I|21|l . A renor- 
malized low energy term, A ren , can be defined by 



A = AnP~ J ^ 



'-du 



(30) 



The free energy is again determined by Eq. 1)25(1 . though 
cannot be evaluated analytically, anymore. The scaling 
behavior of the renormalized tunneling given in Ea. l|29l) 
is no longer a power law, as in the Ohmic case. Still, we 
can distinguish two limits: 

i) The renormalization of A(A) is slow. In this case, 
the integral in Eq. Q25[) is dominated by the region A ~ 
Ao, where the function in the integrand goes as A~ 2 . 
The integral is dominated by its high cutoff, A , and the 
contribution from the region near the lower cutoff, A ren , 
can be neglected. Then, we obtain that -F(Ao) ~ Aq/Ao. 

ii) The renormalization of A(A) is fast. In this case, 
the contribution to the integral in Eq. I|25(l from the 
region A sj Ao is small. The value of the integral is 
dominated by the region near A ~ A rcn . As A rcn is 
the only quantity with dimensions of energy needed to 
describe the properties of the system in this range, we 
expect that -F(Ao) « A ron . 

In the scaling limit, Ao/Ao <§C 1, the vàlues of the 
two terms, A rcn and Aq/A , become very different. In 
addition, there are no other energy seales which can qual- 
itatively modify the properties of the system. We thus 
conclude that only the two terms mentioned above will 
contribute to the free energy. Hence, we can write: 



F(A ) ~ max A r 



A 2 
Ao 



(31) 



The above equation is now used to discuss the possible 
transition between underdamped to overdamped oscilla- 
tions for non-Ohmic environments. Notice that it also 
applics for Ohmic baths. 



Súper- Ohmic 



In the super-Ohmic case (s > 1), Eq. (|30() always has 
a solution and, moreover, we can also set the lower limit 
of the integral to zero. This yields 



A rcn = A e" 



■li 



A j( to 



- díü 



A, 



oe 



-a/(s-l) 



(32) 



For a»l we have A rcn <C A , but there is no transition 
from localized to delocalized behavior. 
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Using Eq. (|3l·|l in the super-Ohmic case s > 1, we can 
approximately write: 



-a/(s-l) 

' A 



(33) 



We thus find a transition from underdamped to over- 
dampcd oscillations at some critical coupling strength 
a~ (s - 1) log(A /A ). 

It is finally interesting to note that the scaling analysis 
discussed in Ref. [30( is equivalent to the scheme used 
hcrc. 



2. Sub-Ohmic dissipation 

In the sub-Ohmic case (s < 1), it is not guaranteed 
that Eq. (|3Ü|I has a solution. In general, a solution only 
exists when Ao/Ao is not much smaller than 1. 

The existence of a phase transition in case of a sub- 
Ohmic bath was first proved in Ref. 31]. Whereas the 
relation in Eq. I|3U|I and a similar analysis based on flow 
equations for Hamiltonians |3^ yields a discontinuous 
transition between the localized and delocalized regimes, 
detailcd numcrical calculations suggest that the transi- 
tion is continuous |33| . 

Since there is a phase transition from localized to non- 
localized behavior, there might also be a transition be- 
tween overdamped to underdamped oscillation. In Ref. 
[3^| , this transition was discussed on the basis of spectral 
functions analogous to the discussion of Ref. [35L l36| for 
Ohmic dissipation. It was found that for s > 0.5 the tran- 
sition takes place for lower vàlues of a as in the Ohmic 
case, e.g., for s = 0.8 and Ao/Ao = 10 the transition 
coupling strength is a* ta 0.2. 

Using Eq. (J301) and Eq. (|2U yields for the sub-Ohmic 



A o /A o 



DL(coh.) 



/ DL(incoh.) 




a 1 

FIG. 2: (Color online). Schematic picture of the difïerent 
regimes in the sub-Ohmic dissipative TLS studied in the text. 
DL stands for delocalized phase, while L denotes a localized 
phase. The lower blue line denotes the continuous transition 
studied in Ref. Is?! . The red line marks the boundaries of a 
regime characterized by a small renormalization of the tun- 
neling rate, Eq. 1301 . and coherent oscillations. 



The scaling of (a x ) is 
d(a x 



dA 



(37) 



The fact that the scheme assumes a fully coherent state 
as a starting point implies that A is not renormalized. 
Inserting Eq. into Eq. l{57|). we find: 



d(<7 x 
dA 



A 



Ao 
A 2 



(38) 




delocalized regime A 2 - ~ 1 
localized regime ^ -C 1 



(34) 



The analysis used in the previous cases leads us to expect 
coherent oscillations in the delocalized regime. 

We can extend the study of the sub-Ohmic case to 
the vicinity of the second order transition described in 
Ref. |3jj , which in our notation takes place for a = 
sAo/Ao <C 1. In this regime, which cannot be studied us- 
ing the Franck-Condon like renormalization in Eq. (|30ll . 
we use the renormalization scheme around the fully co- 
herent state proposed in Ref. |3Ï} ■ To one-loop order, the 
beta-function for the dimensionless quantity (expressed 
in our notation) k = (oA) / A then reads 



-SK + K 



-2 



(35) 



Near the transition, in the delocalized phase, k thus 
scales towards zero as 



k(A) = k [ — 



(36) 



If we calculate (a x ) from this equation, we find that 
the resulting integral diverges as A — > for s < 1. This 
result implies that (a x ) <C 1. For sufficiently low vàl- 
ues of the effective cutoff, A, the value of (a x ) can be 
calculated using a perturbation expansion on Ao, lead- 
ing to (<r x ) ~ A /A . This result implies the absence 
of coherent oscillations, as in the similar cases discussed 
previously. 

A schematic picture of the regimes studied for the sub- 
Ohmic TLS is shown in Fig. 0. 



IV. SUMMARY 

In this article, the entanglement properties of dissipa- 
tive systems were investigated on the basis of the von 
Neumann entropy. 

We first investigated two integrable dissipative quan- 
tum systems -the free dissipative particle and the dissi- 
pative harmònic oscillator - and calculated the von Neu- 
mann entropy. In the former case, this could be done 
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exactly and no non-analyticity was found. The case of 
the harmònic oscillator is the more interesting one since it 
exhibits a transition from underdamped to overdamped 
oscillations for increasing dissipation. This transition is 
also manifested in the entropy, or equivalcntly in the en- 
tanglcmcnt which was calculated in the limit of large bath 
cutoff. 

We also calculated the von Neumann entropy for the 
spin-boson model on the basis of a scaling equation for 
the free energy. Only in the Ohmic case, the resulting 
integral could be evaluated and we analyzed the non- 
analyticity at the transition from underdamped to over- 
damped oscillations. We found that the non-analyticity 
more pronounced than at the actual phase transition. 

In the non-Ohmic case, we argued that the transition 
between coherent and decoherent oscillation takes place 
when the value of (a x ) becomes comparable to the result 
obtained using a perturbation expansion in the tunnel- 
ing matrix, A (as is the case for Ohmic dissipation). In 
the super-Ohmic case, this always yields a critical cou- 
pling strength at zero temperature which differs from the 
analysis in Ref. (3^ . 

In the sub-Ohmic case, the scaling approach can only 
be trusted when the tunnel matrix element is of the same 
order of magnitude as the cutoff. Then a transition be- 
tween coherent to non-coherent oscillations is possible be- 
fore the system becomes localized. For the regime where 
the cutoff represents the largest energy scale, we applied a 
novel renormalization scheme proposed in Ref. |37| |. We 
find that, in the delocalized phase, the system is most 
likely incoherent. 

Concerning the entanglement properties for the non- 
Ohmic case, we were not able to discuss possible non- 
analyticities since the regime is analytically not accessi- 
ble. Numerical work in this direction is planned for the 
future. 

To conclude, we suppose that entanglement properties 
are closely connected to the transition of coherent to in- 
coherent tunneling. Our observations might be uscful for 
future quantum bit manipulations. 
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APPENDIX A: CALCULATION OF THE FREE 
ENERGY OF THE DISSIPATIVE TLS 

We calculate the free energy of the dissipative two level 
system following the scaling approach discussed for the 
Kondo problcm in Refs. \1Í |39i . and formulated in a 
more general way in Ref. |4Çj. For the general long- 
ranged Ising model, the scaling approach was first applied 
by Kosterlitz [IçJ. 



Ti ^2 ^3 (3 

FIG. 3: (Color online). Sketch of the instanton pairs which 
renormalizes the calculation of the free energy of the dissipa- 
tive TLS. 



The partition function of the model can be expanded 
in powers of A 2 as 




(Al) 

where /[(rj — 7j)/t c ] denotes the interaction between the 
kinks located at positions Tj and Tj . A term in the series 
is schematically depicted in Fig. [3]. The scaling pro- 
cedure lowers the short time cutoff of the theory from 
r c to t c — dr c . This process removes from each term 
in the sum in Eq. IjAlfl details at times shorter than 
t c — dr c . The rescaling r c — > r c — dr c implies the change 
A — * A(I + dT c /r c ). The dependence of /[(t^ — tj)/t c ] 
leads to another rescaling, which can be included in a 
global renormalization of A [Ï2l I3H lioj . In addition, 
configurations with an instanton-antiinstanton pair at 
distances between r c and r c — dr c have to be replaced 
by configurations where this pair is absent, as schemat- 
ically shown in Fig. The number of removed pairs 
is proportional to dT c /r c . The center of the pair can be 
anywhere in the interval < r < j3. The final effect is 
the rescaling: 

Z -> Z (1 + A 2 /3dr c ) (A2) 

Writing Z as Z = e~@ F , where F is the free energy, Eq. 
(|A2|) can be written as: 

r)F 

-^ = A*(T C ) (A3) 
In the Ohmic case, the dependence of A on t c = A _1 is 
A(A) = A„(A)° (A4) 
and, finally, we find the following relation: 



dF 


[A(A)- 


2 








A 






VAo) 



(A5) 



This equation ceases to be vàlid for A ~ A ren . For finite 
temperatures, we obtain 

F(T) = j T ^dA. (A6) 
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It is interesting to apply this analysis to a free two level and, finally: 
system. The value of A does not change under scaling. 
We find the following expression: 

^((^) 2 Ao«A (A7) ^) = Z = {tï°fZ ( AQ ) 



9A^\ò iW A«A [A " XX/ dA \1 T«A 

Inserting this expression into Eq. (|A6(I . we obtain 

{ A 2 in qualitative agreement with the exact result (a x ) = 

f ^°<V (A8) tanh(A /T). 
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